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ABSTRACT 


If  the  pairs  of  molecules  in  a  dilute  gas  interact  -with  an  intermolecular 
energy  of  the  form,  CD  =  d/r^  ,  the  collision  integrals  which  determine  the 
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transport  properties  can  be  written 
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Here  u  is  the  reduced  mass  of  the  colliding  molecules.  The  following  constants 
are  evaluated: 


Repulsive  Interaction  Potential,  d  positive 


A(1>(2)  =  0.397601 

A(Z)(2)  =  0.527843 

A(1)(3)  =  0.3115 

A(Z)(3)  =  0.3533 

Attractive  Interaction  Potential,  d  negative 


A(1)(2)  =  0.806907 

A(Z)(2)  =  0.710970 

A(1)(3)  =  0.6750 

A(3)(3)  =  0.4641 

The  constants  for  „  =  2  were  readily  evaluated  in  closed  form.  The  constants 
for  the  repulsive  potential  with  f  =  3  were  evaluated  using  an  expansion  develop¬ 
ment  of  Mott-Smith.  The  constants  for  the  attractive  potential  with  $  =  3  were 
evaluated  numerically  making  use  of  elliptic  integrals  for  the  angles  of  deflection. 
For  attractive  potentials  a  singular  point  was  encountered  corresponding  to  orbiting. 


This  work  is  being  done  under  the  Office  of  Naval  Research  Contract  N7-onr-28511, 


J  O.  Hirschfelder,  C.  F.  Curtiss,  and  R.  B.  Bird,  "Molecular  Theory  of  Gases 
and  Liquids"  (John  Wiley,  1954),  p.  547. 
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The  transport  coefficients  for  a  dilute  gas  made  up  of  molecules  which  have 
a  potential  energy  of  interaction  of  the  form 


( IT)  ~  eS  r 


-% 


(1) 


can  be  expressed^  in  terms  of  the  constants  A  (  F)  ,  when  d  is  positive  and 

^  (J)) ,  \ 

the  molecules  repel  each  other;  or  in  terms  of  the  analogous  constants  (d)  > 

to  be  defined  belowj  when  d  is  negative  and  the  molecules  attract  each  other.  In 
this  paper,  the  values  of  A(l)(%)  ,  A^C?-)  ,  f  S')  and  C^Vs") 

calculated  for  K  =  2  and  S'  =  3. 


are 


Although  the  second  virial  coefficient  would  be  infinite  in  each  of  these  cases, 

the  transport  coefficients  are  finite.  The  evaluation  of  the  collision  integrals  is  not 

difficult  for  either  of  the  repulsive  cases.  When  S"  =  2,  A0'  ( Z)  and  A ^  (2) 

3 

may  be  expressed  in  closed  form.  When  S"  =  3,  a  formulation  of  Mott-Smith  may 
be  used  to  obtain  a  rapidly  converging  infinite  series  for  these  collision  integrals. 
The  integrals  resulting  from  the  attractive  potentials  are  somewhat  harder  to  evalu¬ 
ate  because  of  the  spiraling  of  the  molecules  in  collisions  where  the  impact  para¬ 
meter  is  less  than  a  critical  value. 


1 

2 


J.  O.  Hirschfelder,  C.  F.  Curtiss,  and  R.  B.  Bird;  "Molecular  Theory  of  Gases 
and  Liquids";  (John  Wiley,  1954);  p.  546. 
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We  find  that  for  either  attractive  or  repulsive  point  forces,  the  A  ( I )  and 
are  infinite.  The  second  virial  coefficients  (MTGL,  Ref.,  1,  p.  156)  are  infinite 
for  either  attractive  or  repulsive  potentials  if  £  is  not  greater  than  3. 
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H.  M.  Mott-Smith,  "A  New  Approach  in  the  Kinetic  Theory  of  Gases",  M.  I.  T. 
Lincoln  Laboratory,  Group  Report  V-2  (December,  1954). 
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When  %  =  2,  this  spiraling  results  in  an  infinite  angle  of  deflection  whenever 
the  impact  parameter  is  less  than  the  critical  value,  and  causes  singular  behavior 
in  a  part  of  the  collision  integral.  This  part  of  the  integral  is  approximated,  and 
the  remainder  is  obtained  in  closed  form. 

For  'b  =  3,  the  critical  impact  parameter  corresponds  to  a  singular  point 
in  the  collision  integrals.  The  asymptotic  behavior  in  the  vicinity  of  the  singular 
point  is  investigated  and  the  integrals  are  evaluated  numerically.  Use  is  made  of 
the  fact  that  with  5=3,  the  angle  of  deflection  is  simply  ekpressable  in  terms  of 
incomplete  elliptic  integrals  of  the  first  kind. 

4 

The  collision  integrals  are  defined  in  the  following  way.  The  collision  of 
two  molecules  is  considered  as  the  equivalent  one  body  problem.  A  particle  with 
reduced  mass  n  =  approaches  a  fixed  center  of  force  with  initial 

"  rry\ ,  4 

speed  g  and  impact  parameter  b.  The  potential  energy  of  interaction  is  given  by 
Eq.  (1),  where  d  is  a  constant,  and  r  the  distance  of  the  incoming  particle  from 
the  scattering  center.  Figure  1  shows  the  trajectory.  The  particle  is  deflected 
through  an  angle  of  deflection  ^  and  the  distance  of  closest  approach  is  r 

m 

Now,  if  we  define 


the  angle  of  deflection  may  be  written  as  a  function  of  y 


Xfyj  = 


-  rr 


=  IT  —  2  ^  [  I  -  %Z~i(  ]  ciy  (3) 


MTGL,  Ref.  1,  p.  546.  It  should  be  noted  that  our  A  a)  and  that  of  MTGL 
are  equal  to  2  ^  A/s  where  the  /\j(5+l)  are  the  collision  integrals  of 

S,  Chapman  and  T.  Cowling,  "The  Mathematical  Theory  of  Non  Uniform  Gases" 
(Cambridge  Press,  1939),  p.  172. 
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Fig.  1  Trajectory  for  repulsive  potential. 


Fig.  2  Trajectory  for  attractive  potential. 

Case  where  particle  passes  through  center  of  force. 


-4- 


The  upper  limit  on  the  integral  is  the  first  positive  root  of  the  polynomial  in  the 
denominator  of  the  integrand. 

Then  the  collision  cross-sections  for  the  various  transport  properties  are 
given  by 

/V"  _  ^  <r  r  jjl  rd  Atrjtz) 

(4) 


Qtri  =  A't>is) 


tj  \  _ 

where  the  molecular  collision  integral,  A  ( £  )  is  given  by 

Amd)=  /“r  I-  c^>x)y.  Jy* 


(5) 


When  the  potential  is  attractive,  d  =  -a  and  the  potential  has  the  form 

(6) 


=  -  o- 

Also,  let  us  define 


(7) 


It  is  easily  seen  that  with  A  ,  which  is  a  real,  positive  quantity,  Eq.  (3)  becomes, 
for  an  attractive  potential 


?**(/*> 


/y 

Xy'rA'-P'  r  -j  -/j 

[  (8) 

If  ^6  is  greater  than  0  cr^i  the  polynomial  in  the  denominator  of  the  integrand  has 
a  first  positive  root,  ym(/3  )•  However,  if  {3  is  less  than  0  ,  y^  (^3  )  =  CO 

corresponding  to  the  fact  that  the  incoming  particle  actually  passes  through  the 
center  of  force,  r^  =  0,  and  then  goes  out  again.  In  this  case,  the  angle  of  deflec¬ 
tion  is  given  by 

r  °°  y 

■x(fi)=  -29,  =  -2  f  [,-f  +  Hy)  ]  J 


#  (9) 


rather  than  by  eq.  (8).  Such  a  trajectory  is  shown  in  Fig.  2. 

For  the  attractive  potentials  the  various  collision  cross-sections  are  given 
by 


(10) 


i 


where 


cwii)  =  S,  ( 1  - 


I.  Centers  of  Repulsion 

For  repulsive  potentials,  when  £  is  2  or  greater,  the  A  \  are  finite. 

The  case  of  S  =  1,  however,  is  a  potential  of  sufficiently  long  range  that  the  integral 

(_p) 

foh  A  ( /)  is  infinite.  In  this  case,  the  angle  of  deflection  is 


X  (y.)  =  2  ^uc'  [(t 


and  the  integral  in  Eq.  (5)  diverges. 

(a)  F  =  2. 

The  repulsive,  inverse  square  potential  is  of  interest  because  the  integrations 

5 

may  be  carried  out  analytically.  First  of  all,  the  angle  of  the  deflection  from 
Eq.  (3)  is 

/y~  _  y 

[  1  ~  (  1  *  Jf  (13) 


Let  x  =  2  yt  ;  then  y^  =  (1  +  £  )”  'i 


%  (x)  -  rr  -  2  ( i  +  i )~K 


(  l+  x) 


x  Cx)  =  rr  [  i  -  i 


The  evaluation  of  A^  ^(2)  and  A^  ^(2)  was  carried  out  by  John  S.  Dahler  with  the 
help  of  Donald  W.  Jepsen. 
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With  this  change  of  variables,  Eq.  (5)  becomes 

2)  ~  f  (  I  ~  Jx 


(15) 


Z  2 

Now  let  x  =  z  /  (  1  -  z  )  ,  then  %  =  TT  (  1  -  z).  Substituting  into  Eq.  (15) 
and  integrating  by  parts 


(1) 


A 72)  =  -t 


_  i—  -f 

“  O  I 


%  /' 


77*  Z 

I  -  z1 


Jz 


(16) 


-£  r'  jz 


I  ~  Z‘ 


(17) 


But  it  is  easy  to  show  that 
, 1 


^  z  Jz  —  pf  f-  0  ^  5 1  (0.  ^  V')  —  2  2c  (<r\ir)-^ 


(18) 


where 


=  f 
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(19) 


Thus  we  obtain 


A(1)  (2)  =  0.  397601 


A*2)  (2)  =  0.  527843 


(20) 

(21) 


"Tables  of  Sine,  Cosine,  and  Exponential  Integrals",  Vol.  II,  WPA  (Federal 
Works  Agency,  Works  Progress  Administration  for  the  City  of  New  York. ) 
(1940). 
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(b) 


=  3. 


In  the  'case  of  the  inverse  cube  repulsive  potential,  we  were  able  to  make  use 

3  <£) ,  0 \ 

of-  a  method  which  has  been  developed  by  Mott-Smith  to  evaluate  the  A  (  o  ) 

7 

Through  the  application  of  Btlrmann's  Theorm  he  has  found  that  Eq.  (5)  can  be 
written  as 


A'*(V 


rir 

=  5TyS  J  (l-e<r*VK)  F(x)^y i  ^  J% 


(22) 


where  the  F  (  X  )  sinX  is  given  by 


co 

—  (bcrt  l'X  ft’ 


->  -  Vr 


/°=0 


The  are  constants,  depending  only  on  5"  • 


(23) 


D, .=  (24> 

’  '  K  Jj!=6 


with  the  H(  2  )  given  by 


/-/( 2)  " 


oo 

z 

5-/ 


“A 


t:0 


+j)T(-s  +  i) 

r  a  st  -  s  *■!> 


(25) 


For  =  3,  we  obtain  the  results  given  in  Table  I. 


Table  I 

The  Mott-Smith  Coefficients  D^>  for  f  =  3 


p 

D/° 

0 

0.  25438 

1 

0. 12487 

2 

0.01019 

3 

-0.00282 

4 

0. 00317 
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E.  T.  Whittaker  and  G.  N.  Watson:  "A  Course  of  Modern  Analysis";  4th  Ed, 
(Cambridge  University  Press)  (1952)  p.  128. 
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Substitution  of  Eq.  (2.3)  into  Eq.  (22),  and  integration  gives  the  general 
expression 


/fYr)  - 


^  ! 

)_ 

p~0 


(26) 


The  two  cases  which  are  of  interest  are  then 

co  i 

An)(3 )  -  ft  £  P  +  %]. 

Xi](3)  =  Dp[A% 

p~e> 

Thus  we  obtain 

A(1)  (3)  =  0.  3115 
A(2>  (3)  =  0.  3533 


(27) 


(28) 


(29) 

(30) 


II.  Centers  of  Attraction 

When  the  potential  is  attractive  there  is  a  singularity  in  the  expression  for  the 
angle  of  deflection.  The  Mott-Smith  formulation  is  no  longer  applicable  so  we 
proceed  as  follows. 

There  exists  a  critical  value  of  b  =  bc  for  any  relative  kinetic  energy 
S  ^2  such  that  for  b  less  than  b^^  there  are  no  positive  roots  to  the 
equation 

f  (y)  =  /  ~  ft  +  “  0  (31) 

where 

I  __  _ 


c 


(32) 
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and  for  b  =  b  .  ,c  =  c  ..  This  critical  value  may  be  obtained  by  noting 
crit  crit 

that  Eq.  (31)  has  a  double  root  at  the  critical  value.  Therefore,  both  f  { y )  and 
f'  (y)  vanish  for  the  same  positive  root,  y  * 


f V-y)  =  r‘ +  c  ~  ° 


(33) 


One  may  combine  Eqs.  (31)  and  (33),  to  solve  for  this  root 
^  6  ' 


f  fk]'1 


and  then  go  back  and  obtain  the  critical  value  of  c 


^  CW-vt 


2  S'H  [  S-2]k~' 


Thus 


-‘/t 


=  2  l  J-a 


(34) 


(35) 


(36) 


(a)  £  =  2. 


For  the  attractive  case,  with  §  -  2,  the  angle  of  deflection  is,  from 


(37) 


and 
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Thus, 


X  ( fit  -  IT  [  I  -  (I  ip  )  3  ]  ;  /3  '  >  1 

a  (p)  =  -2  (& 


2-  ,  / 


In  spite  of  the  infinite  value  of  the  angle  of  deflection  for  all  values  of  {£ 
less  than  ,  one  may  still  obtain  the  collision  integrals.  Eq.  (11)  is  split  into 
two  parts 


-(a) 

C  w 


J  (  j  -  OrtMz)/2  +  J  k  (l- 


-  G  +  H  (4D 

The  second  of  these,  the  H  ,  are  obtained  in  closed  form.  Taking  the  expres¬ 
sion  for  from  Eq.  (39),  and  letting  x  =  (  1  -  )  ,  it  is  found  that 


f  oo*sjF*  c 0(^hr/) 


!  jrCO  f\  f  / 

y  J  ju/a  2  Trx  &  (  ^4-/ 


These  expressions  may  be  integrated  by  parts  to  give 


ji>sv\  fry 

l  -  xz 


t  r  ^  ^ 


(45) 
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These  may  be  rewritten  as 


H' 


f  C 03  yi^sA  fr  X 

Ay  — 

C 1  TP  X 

l  Jo  !  -  X* 

CX  r- 

y- 

i 

X 

> 

(46) 


„(2)  _  If 

H  -  ^ 


-J.  -T^r  +  i  -T=ir*  ‘Vxj  <«>' 


The  first  integral  in  each  case  is  found  in  Bierens  de  Haan  ,  and  second  we  found 
in  the  evaluation  of  A^  (2)  and  A^  (2).  Thus 

f  ^  **» 

H(1)  =  f :j  [Ct(rr)^J«7r  -  St  (ft)  wTr]  *2  m 


H 


(2) 


=  -f( -[Ci(w)^a*  -Si(*ft-)c<n*tr]  +i  (49) 


Therefore, 


h(1)  =  SL  (zrr) 


H<2)  =  Si  Otn) 


0.  556907 


0, 585970 


(50) 


(51) 


(  Jt  ) 

The  values  of  the  integrals  G  may  be  approximated  in  the  following  way. 


AJ  ) 


j  f  7  (  I  -  c<xJ%)  J  p1 


(52) 


D.  Bierens  de  Haan,  "Nouvelles  Table's  D'Integrales  Definies",  (G.  E.  Stechert 
and  Co.,  1939),  Table  1 6 1 ,  3. 

J  frrj*  J/  ~  i  f  C Hpi)  ^  PI  " 

*  ft  r  V  J 


With 


rP 

Cl  (p$)  = 


PI  X 


Jx 
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If  the  upper  limit  on  the  integral  for  the  angle  of  deflection  is  taken  to  be  an  arbi¬ 
trarily  large,  but  finite  number  2> 


X 


%  Up-  ~i)  tl*  ~ < 


(53) 


With  this  expression  for  the  angle  of  deflection,  the  integrand  of  Eq.  (52)  is  a 
very  rapidly  oscillating  function  of  (2^1  oscillating  between  0  and  2  for  , 

and  between  0  and  1  for  (j  .  As  ~2L  becomes  very  large,  the  peaks  become 
very  close,  and  the  average  value  of  the  integrand  becomes  very  close  to  1  for  the 
and  to  Jp  for  the  •  Thus,  in  the  limit  that  2  becomes  infinite,  cor¬ 

responding  to  the  actual  kinetic  theory  problem, 


G";  = 

0,2  5* 

(54) 

g(3)  = 

0 . 1  25" 

(55) 

Thus  we  obtain 

C  (*)  = 

G,"'  *  H0> 

=  0.  806907 

(56) 

C  uVs)= 

Gu>  +  Hll) 

=  0.710970 

(57) 

(b)  %  =  3 


In  the  case  of  the  inverse  cube  attractive  potential,  ^ 


crit. 
'A 


)1/6  and 


y  (/3  )  varies  from  unity  for  P,  -  oo  to  y  (  &  )  =  3  *  .  Because  of  the 

Jm  r  r  7m  I  crit.  .j  . 

singular  point  it  is  convenient  to  divide  the  integration  of  C  (3)  into  two  parts 


C  (3)  ~ 


where 
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(£>) 

C,  '*> 


1 

2 


-(?) 


^3 


(59) 

(60) 


^  [  I  -  c^L  (20,)]  J p7- 

C"'(  3)  =  i  Li  1  ~ 

(?) 

Here  Q  and  Q  are  given  in  Eqs.  (9)  and  (8)  respectively. 

X  Lt 

( a ) 

In  evaluating  C  (3)  it  is  convenient  to  change  the  variable  of  integration 

2  2  1  2 

from  A  to  q  to  q  where  q  is  minus  the  (only)  real  root  of  the  polynomial 
in  the  denominator  of  Eq.  (9) 

ft1  =  f'tlC-f)]  3  (61) 


2  1/3 

Since  the  real  root  is  negative  in  the  range  0  (  f3  ^  (^l)  ,  q  is  positive. 

Making  the  change  of  variables, 


«-  (4>“  P- 


T 


I  -  C, 


% 


P  +  yj  Jy  <62> 


9 

which  is  an  incomplete  elliptic  integral  of  the  first  kind  . 
Thus, 


0,  - 


_lf  A 

I  -  Cl 7 


where 


§•  =  [ 


xx, ->i 


(63) 


(64) 


(65) 


(66) 
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P.  F.  Byrd  and  M.  D.  Friedman,  "Handbook  of  Elliptic  Integrals  for  Engineers 
and  Physicists",  (Springer,  1954),  p.  88,  Eq.  241.00. 
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&{&')  ~S/3  (3-f)  ^(J»  , 

Since  — ■  -  r  —  3  - - — ^—7-  ,  the  first  integral,  f  (3)  if 

j(f>  (hf)* 

now  given  by 


^  (J)  . 

C,  «) 


1 

~//3  I  [  /  - 


n-%1). 

a-fy 


Figs.  3  and  4  show  the  integrands  for  J)  -  1  and  for  J  =  2  respectively 
as  functions  of  q^. 

Z 

It  is  seen  that  the  integrands  oscillate  wildly  as  q  approaches  0.  75.  The 

2  2 

integrals  are  subdivided  into  four  ranges  of  q  .  For  the  smaller  values  of  q 

numerical  integration  is  satisfactory  making  use  of  the  elliptical  functions.  How- 

2 

ever,  for  the  larger  values  of  q  it  is  convenient  to  use  2  (7,  as  the  integration 

J  f  ft* ) 

variable  and  develop  approximate  expressions  for  which  will  allow 

analytic  evaluation  of  the  required  integrals.  Thus,  C^^^(3)  is  divided  into  the 
ranges  0  -  q^  ~  .  52,  .  52  ^  q^  .  68,  .  68  $  q^  K  .  748,  and 

.  748  3  q  ^  .75. 


^  (j)  T  (J)  ,  -T-U)  -r-fj)  ,  -1-W) 

C,  73)  --  Z  +  J,  +  I,  +  T, 


where 


l!"  = 


[l- 


J(sa,)l  -3 ~  $ 


'S,  91  2U 


£  I  —  (  2.  0  1  )  J 


M£) 

JU&) 


J(29t) 


,h\  rll>,Vt*l(s  At  A*) 

l3  =  i  friol„  r;  -  jfa}  j(w)  aD 

jwl ,  j  r  i  1-  im  (72) 

H  *  JU  ,  VU6>  *  ' 


.15 

3.0  .25  .30  35  MO 

.45  .50  ,55  .60  .65  70  .75 

Fig. 

3  Integrand  of  C  ^(3) 

from  Eq.  (67) . 

0  .05  .10  ,15  .20  .25  ,30  .35  .40  .45  .50  55  .60  .  65  .70  75 


(2) 

Fig.  4  Integrand  of  '(3)  from  Eq.  (67). 
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The  integrals  1^  and  were  evaluated  numerically  using  Simpson's 

2  2 

Rule  with  interval  sizes  of  .04  in  q  from  0  £  q  £  .  40  and  intervals  of  .02 

2  2 
for  larger  values  of  q  .  The  integrands  were  evaluated  at  these  values  of  q  by 

making  use  of  Eqs.  (63)  and  (67)  together  with  interpolated  values  of  the  elliptical 

integral.  Four  point  Lagrangian  interpolation  coefficients  were  used  to  interpolate 

along  both  the  rows  and  the  columns  in  the  table  of  Byrd  and  Friedman.  ^  Thus  we 

obtained 


(1) 


.  2449 


(73) 


(2) 


.  1136 


(74) 


An  expression  for 


MF) 


d  (2  &i) 

ained  by  fitting  the  cm 

,  with  a  cubic  equation.  (See  Fig.  5  for  a  plot  of  ^ 


,  sufficiently  accurate  in  the  interval 

2  ^2 


5.87292  £  ZQ1  .£  9.  10288  is  obtained  by  fitting  the  curve  of  .  —  f3  vs. 


crit. 


crit. 

-  ^  2  vs. 


2 


Pi*  -  P  =  D  +  tz  (3.6,)  4  F  (26,?  4 Q  (26,  f  (75) 


The  points 


201 

5. 87292 

6. 74996 

7.94016 

9.  10288 

A  2  2 

P  crit  P 

.  500778 

.411380 

. 300569 

. 209800 

2 

q 

.  52 

.  58 

.64 

.  68 

are  used  to  evaluate  the  constants  in  Eq.  (75).  The  constants  are  D  =  1.062129, 
E  =  -.0640403,  F  =  -.0092255,  and  G  =  .0006563.  With  these  constants, 

Eq.  (75)  fits  the  actual  curve  better  than  .05%  in  the  region  of  interest. 


10 


Ref .  9,  pp.  324. 
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Integration  of  Eq.  (70)  gives 


I 


(i) 


.-r(l) 


_  ± 
2. 


We  obtained 


Efa0,)  +  Fpft)1  +  G 

~  fj;  4-  5  F  (30,)  4  3  C?  (*?$>)'  ■  Q  (~f  ~J  ju-h l 

-[if  +  6G  (Z&,))  otn.  be,) 


2ft  =  V03LW 


-£{  E(2B.)  +  F (29,)*  +  G  (26,? 

4  7j  [  E  +  2F(M)  +3G-(29,)3-}G}jj^m,) 
+  i{  F  +  3GM)J  C*(V6,) 


(76) 

X9,=  5“,  f  7292 

x&t~  9,  /ca?# 


J. 


(77) 


4jg),  -  F,  8-72  71 


I2  (1)  =  .1020  (78) 

I2  (2)  =  -0690  (79) 

tjj )  (J) 

The  integrals  and  I  are  evaluated  with  the  aid  of  an  asymptotic 

expression  for  d  (p)  f  d(2& accurate  near  cr^  •  ^he  asymptotic  form 
for  Q  ^  as  approaches  .  75  may  be  obtained  as  follows: 


From  Eqs.  (63)  and  (64), 

B,  =  ?f o-tVa-f)]'*  F  (  P,k) 

(80) 

*\  i  -l 

As  q  ->.75,  k— >1  and  120°  .  Since 

F(¥>,k)  =  2K(le)  -  T  (rr-y’tk)  , 

(81) 

VI f,k)  cs  lK(k)  -  T((>o°,  I )  . 

(82) 

11 


Ref.  9,  p.  12,  Eq.  113.02 


/  /  2Y* 

If  k  =  [  I  -  k  I,  the  expansion  for  K  (k)  is 


K(k)  =  (J*$)  [  I  +  ¥  k'‘]  -  ?  >' 


-jX  I  t  i 


Since  k  —  1, 


2  K  -  2  Jn  k‘  -  YY  /-*“ 

r  ( v,  k)  £  Jm  jrfc  -PUo'O 

2  -a  2 

If  the  substitution  q  =  -  z  is  made, 

T 

-  /-  f  ^  +  ■" 


/-k" 


i4  [  /  *  Vz*  +  *  •  '  J 


Therefore, 


e,=  jl  it 


-  T(to%k) 


Since 


-TUo\\)--  -Jm  (hv  6a*  U')*J»(2- i*)  m 

2*  =  3  i  (-  ~ 3  ‘)  Q  .  *  (90 


Ref.  9,  p.  298,  Eq.  900.05 


See  Appendix  A  for  an  alternate  derivation  of  this  equation. 
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~~ 

Jtfi1) 

J(z') 

jua)  ~ 

J  (zx) 

d  ( <2  &■/ ) 

(91) 

From  Eq.  ( 6 1 ) ,  it  is  found  that  asymptotically 

4 

[  1  -  £  Z: 

'•] 

(92) 

Therefore, 

MA) 

<*(!&.) 

=  3C(^ )*  (s- 

17,  e 

-  k  c~ 6' 

(93) 

_ 

tS  a 

J+~ 

II 

p  oo 

2k  f  e 

~9,  J  &, 

(94) 

/C  -  f  - 

-  xl  e~ 

e, 

(95) 

This  suggests 

approximating  ^ 

2  , 

by  an  equation  of  the 

form 

i 

p  crit. 

I 

r° 

n 

i 

e'[A  +E 

e  &<$,)'] 

(96) 

over  the  range  9.  10288  ^  2  6^  ~  16.91616.  The  points 


2^1 

9. 10288 

11. 74086 

16.91616 

/S'crit^2 

. 209800 

.0780498 

.0072236 

2 

q 

.68 

.  727 

.  748 

are  used  to  evaluate  the  constants  in  Eq.  (96).  The  constants  afe  A  =  -30.419284, 
B  =  7.524143,  and  C  =  -.219515.  Eq.  (96)  fits  the  actual  curve  better  than  .5% 
in  the  region  of  interest.  Eq.  (71)  may  now  be  integrated  analytically;  the  results 


•  21- 


are: 

_  ( i> 

— L  a 


_  -9, 

e 


O.^f/4  +  8(29,)  +CCX6,f} 

-  {(6,1  A  -0.3 2B  +0,12  SC) 

+  (6 ,IB  -  6,Q  ic)(7  9,)  ■+  6, 1C  lie,)1}  cM®)\ 
1  {  (a.lX  -  CMS  ~0,l0tC) 

+  (0,2H  ~0,UC)(lb')  +  6,2  (2&) 


26,-- U,  VC  16 


(94) 


l9,*<},/0i  6? 


lT=  -e* 


(Q.b^lU  A  -  (),  assziMB  +0.68¥C  79 C) 

-+(b.  WlOiB-OMonzbC)^,)  1  6,%47D6Cm,)L' 

j  +  {(0,0  6~  StJVA  ~0r/03ge}6/3~6,67CS32C) 

+(6,  6sem&  "  0,267  c,  ISC)  (fz)  +  (6.  oft  Me, )(l6>,)Z]^0ij\ 
-[(c.awiA  -6,110717  B  tC,ICj3¥6C) 

+  (0,02Vt/lB-  6,221  *S¥C)  (29, )  +6,677  me  {X6,)’L]  c*xM) 


16,=  !L,  vac 


s6izeUoi$v 


(95) 


Thus  we  obtained 


I3(1)  =  .  1217 


I3(2)  =  .0444 


(96) 

(97) 


(1) 


The  integrals  I4'  '  and  are  e valuated  with  the  help  of  Eq.  (93). 

14  ~  ~k  C~&'[  I  *  0,9  (  (*&,)  ~  g  Cen(3&,))ll  (98) 

10.  71616 


-r  (2)  i  Q  f  .  _  |  £6,  z  GO 

i,  =  -  «e  ‘[i  -  o.inwtt^m)-^^,)**] 


2  0,-16.11616 


(99) 
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Thus  we  obtained 


I*1*  =  .0026 

4 


I }2^  =  .0022 

4 


(100) 

(101) 


From  Eqs.  (68),  (73),  (78),  and  (100)  and  Eqs.  (68),  (74),  (79),  and  (101) 
respectively,  we  obtained 


C^^{3)  =  .4712 


(102) 


C  (Z)(3)  =  .  2292 


(103) 


To  evaluate  the  (3)  ,  as  defined  in  Eq.  (60),  it  is  convenient  to  change 

2  .  2  2 

the  variable  of  integration  from  p  to  p  ,  where  p  is  the  first  positive  root 
of  the  denominator  in  equation  (8)  , 


p-i)] 


v. 


(104) 


Since  y— ( ^  )  =  p,  it  follows  that 


a-  tp'" 


which  is  another  elliptic  integral  of  the  first  kind. 
Thus 


13 


(105) 


(106) 


where 


^  —  2  [  p  (f  (fyx-3)*] 


\ 


(107) 


13 


Ref.  9,  p.  74,  eq.  234.00. 
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I,'2  =  j  + 


(108) 


.  Vi  +  3 

[  3.(3-?') 


(109) 


Eq.  (60)  for  the  C  ^(3)  becomes 

,  /-3 

^(j>)  -L.fr, 


C\ 3)  --  ^  J  f 1  -  Jp 


(3- pz) 


(110) 


Figures  6  and  7  show  the  integrands  for  J)  =  1  and  J  -  2  respectively,  as 
2 

functions  of  p  . 

Here  again,  the  integrations  were  carried  out  in  several  parts,  with  the 

2 

method  chosen  to  fit  the  behavior  of  the  integrand  in  the  various  ranges  of  p  . 
C^(3)  is  divided  into  the  ranges  1  ~  p^  1.025;  1.0  25-^  p^  ^  1.050; 
1.050  $  p  £  2.90;  and  2.  90  $  p  $  3.00.  Thus, 


~(J>)  .  -r-  b>)  *-i-  Ml  .  -r  (J>  T 

C,  ft)  =  J,  4 1  4  J,  +  J„ 


(  HI) 


where 


-ru)  _L 
J,  - 


f  I  -  (-0* tel* (29*)] 


(112) 


-r  I*'  - - 

J2  ~  3'3* 


(113) 


(114) 


jf-fi  -1 

J  v  r  5  ' ' 


(-l>J  C0*-*(2%)}  (pf#,  Jf> 


(115) 
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To  obtain.  ^  ,  write  Eq.  (10  5)  as 

e,= 

2  a-*  2  v  2 

Since  p  —  1,  the  term  (p  -  1)  (fjp)  can  be  neglected,  and 

f/[i -<*>']'%  -  ?/> 

Using  the  fact  that  [l  -  (-1)^  cos^  (Z  0  ^)]  =  (1  -  cos^X-  )  and  Eq.  (112)  , 

j/J) "  J"  (I"  Jf>'L 


(116) 


(117) 


(118) 


Since  p  is  almost  1,  Eq.  (118)  may  be  written  as 


. ( i-o'x)  s*(f-\)’%jp' 


(119) 


From  Eqs.  (8),  (117),  and  a  power  series  expansion  of  cos  X,  one  finds  that 


I  -  anryt  -  ?  (p-i)7'  +<*  (y- if  -■ t-  ..  . 


(120) 


Powers  of  (  p  -  1)  higher  than  three  are  neglected  in  Eq.  (120)  and  oL  is  chosen 
so  that  the  integrand  of  Eq.  (119)  is  equal  to  1.  15406,  the  value  of  the  integrand  of 
Eq.  (110)  when  p^  =  p  ^  =  1.025.  0(  equals  239.49.  From  Eq.  ( 1 19) ? 


(121) 


=  VCh  \  f  V')%  +?(*“■  **■<)(*->)%  &(?,-/)'* 

-I1  (i  -  h  l*'*  -%+%->>*  (122) 
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Thus  we  obtained 


=  0.0015 


(123) 


Ji(2)  =  0.0029 


(124) 


Four  point  Lagrangian  integration  coefficients  are  used  to  get  the  J„ 


J2(1)  =  0.0026 


(125) 


J2(2)  =  0.0081 


(126) 


(J) 

J  is  obtained  by  using  Simpson's  rule  with  interval  sizes  of  0.025  in 

2  ^  2  2 
p  from  1.0  5  to  1.20,  0.05  in  p  from  1.20  to  1.30,  and  .  10  in  p  from 

1.30  to  2.90 


J3(1)  =  .1996 


(127) 


(2)  _ 


(128) 


^  is  found  by  noting  that  an  upper  limit  to  the  value  of  Eq.  (115)  is 
obtained  by  setting  cos  2  &  =  1.  Thus, 


y  Tti)  -4 

J  */  ^  ~  3 


=  0,00026 


(129) 


Since  the  cosine  term  of  Eq.  (115)  has  infinitely  many  oscillations  in  this  region, 
to  a  sufficient  degree  of  accuracy,  the  factor  (1  +  cos  2©2)  can  be  replaced  by 

its  mean  value  of  unity  in  the  integral  for  J  ^  .  Thus,  a  good  estimate  of  J  ^ 

I  (1)  4  4 

is  -r  J  and 

%  max 


.0013 


(130) 
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The  best  estimate  of  J 


(2) 


is  obtained  in  a  similar  fashion.  Here,  the 


maximum  value  of  (1  -  cos  2  Q  )  is  unity  so  that  the  maximum  value  of  J. 

/  (l)  c  2  4 

~2~  ^Max  '  avera8e  value  of  {1  -  cos  ZQ is  one -half,  so  that  our  best 

estimate  of  J  ^  is  Jr  J  ^  ^  or 
4  ^  Max 


(2) 


Eq.  (HI) 


From  Eq. 


2 

J  '  ;  =  .0007 

4 

(  131) 

gives 

C  ^  (3)  =  .2038 

(132) 

C2(2)  (3)  =  .2349 

(133) 

(58),  the  collision  integrals  are  found  to  be 

C(i)  (3)  =  .6750 

(134) 

C(2)  (3)  =  .4641 

(135) 
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SUMMARY 


The  collision  integrals  we  evaluated  are: 


Repulsive  Potential 


A(1)(2)  = 

.  397601 

A(2)(2)  = 

. 527843 

A(1)(3)  = 

.  3115 

A(2>(3)  = 

.  3533 

Attractive  Potential 

( 1) 

CV  '(2)  =  .  806907 

C(Z)(2)  = 

. 710970 

C(1)(3)  = 

.  6750 

C(Z)(3)  = 

.4641 
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APPENDIX  A 

In  this  appendix,  an  alternate  derivation  of  the  asymptotic  equation, 
Eq.  (90),  is  given.  Write 


e,=0„  '%  +  X  (a-d 


and  make  the  substitution  x  =  1  -3^  in  Q ^  and  z  =  3  5 -1  in 

Then 


$„  =  J^'[x'(i-Hx)  +  czd-*)?]  /%Jx 
Q=  f°[zz(i*y^)  fz]3Y'tJz 

I  3-.  u/ 


(A -2) 


(A -3) 


Cl  =  3*  [  3(3*  ~  3fil^]  . 


/3*Zi  (A-4) 


2 

Let  Xq  be  a  small  positive  number  and  take  ^  close  to  {3  Ciit  so  that  c  is 


much  smaller  than  x^  .  Then 


6n  =  (X°  (  X'1  +  1/2  Jx  +  f  X~‘ (!'%*) 

*0  ■ 


c.  [x*  X,  c  _/  (  I  t  3-/J(f  ,  _  %  x()<i  -/ 


(A -5) 


e„  = 


n  (?■  -  3^ 


(A-6) 
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In  order  to  evaluate  Q  ,  ,  take  z  small  and  choose  much  less 

A.  Cd  U 


than  Zq  .  Then, 


l(z  *  +3  c'za  +c'-)l/x  -h?z0+3c*-) 


+  X 


2  c  +  3  cl 
(H%Z.)K  -/ 

(  /  +  +  77, 


e,  =  ^  t 


(A-7) 


(A-8) 


Therefore, 


<3,  - 


(2-3K) 


C 


(A-9) 


2  3  2 

By  using  the  equation  q  =  ^  z  ,,  we  obtain 


^  A»* 

c  = 


4  z 


I  —  2  Z‘ 


a  s M s 


(A  - 10) 


(A-ll) 


Eq.  (90)follows  immediately. 
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APPENDIX  B 

In  this  appendix,  values  for  the  angle  of  deflection  ~X_  are  given  together 
with  a  graph  of  ^  vs.  ^  . 

2 

The  angle  %  is  evaluated  for  values  of  q  listed  in  Table  a  by  using 

Eqs.  (9)  and  (63).  Since  in  this  case  ^  90°  ,  Eq.  (81)  is  also  used. 

Four  point  Lagrangian  interpolation  coefficients  are  used  to  interpolate  along  both 

14 

the  rows  and  columns  in  the  table  for  F(  ^  k)  found  in  Byrd  and  Friedman  for 
2  2  2' 

all  q  ,  For  all  values  of  q  except  q  =  .734,  .741,  and  .748,  the  four 

point  Lagrangian  interpolation  coefficients  are  used  to  interpolate  along  the  column 

2 

of  the  table  for  K  (k)  .  For  the  three  values  of  q  just  given,  it  is  more  accurate 

,  .  15 

to  use  the  equation 

K  =  A  +  k1  +  (A.  -  £)  }t  +  “  •  (b - 1) 

where 

A  ~  V 
k'  =  0-  ft"1 
2 

For  values  of  p  listed  in  Table  b,  Eqs.  (8)  and  (106)  are  used  to  evalu¬ 
ate  X  .  Four  point  Lagrangian  interpolation  coefficients  are  used  as  explained 
above  to  evaluate  the  elliptic  function. 


(B-2) 

(B-3) 


14 


Ref. 


10. 


15 


E.  Jahnke  and  F.  Emde,  "Tables  of  Functions",  (Dover,  1945),  p.  73. 
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Table  a 


2 

crit 


) 


b2 

2 

q 

Bl 

(3  2 

2 

q 

0 

0 

0 

0.40778224 

0,  52 

5. 87292 

0.01976052 

0.04 

1. 14281 

0.43565336 

0.  54 

6. 14246 

0. 04065843 

0.08 

1. 64778 

0.46537752 

0.  56 

6. 43332 

0. 06282202 

0.02 

2.  06022 

0. 4971807 

0.  58 

6. 74996 

0. 08640116 

0.  16 

2.  43106 

0. 5313292 

0.60 

7.  09994 

0.  11157215 

0.20 

2.  78204 

0. 5681396 

0.  62 

7.  49232 

0. 13854409 

0.  24 

3.  12378 

0. 6079917 

0.  64 

7.  94016 

0. 16756714 

0.28 

3.  46410 

0. 6513440 

0.66 

8.  46456 

0. 19894356 

0.  32 

3.  81040 

0. 6987600 

0.68 

9.  10288 

0. 23304244 

0.  36 

4.  16814 

0. 7509361 

0.  70 

9. 93132 

0. 27032007 

0.40 

4.  54440 

0. 8087475 

0.  72 

11. 13542 

0. 29032412 

0.42 

4.  74186 

0. 8305105 

0.  727 

11.  74086 

0. 31134829 

0.44 

4.94660 

0. 8531540 

0.  734 

12.  55009 

0. 33348871 

0.46 

5.  16038 

0. 87674009 

0.  741 

13.  79346 

0. 35685478 

0.48 

5.  38462 

0.90133673 

0.  748 

16. 91616 

0. 38186008 

0.  50 

5. 62188 

0.90856029 

0.  750 

oo 

-33- 


Table  b 


</?2 


/32 

2 

P 

-7< 

£ 

2 

P 

-  X. 

0.. 908561 

3.  000 

- —  -  - . - 

CO 

1. 08128 

1.  600 

1. 22953 

0.908832 

2.  900 

7.  33231 

1. 14471 

1.  500 

1. 01080 

0.909693 

2.  800 

5. 88405 

1.  23977 

1. 400 

0.80593 

0.911285 

2.  700 

5.01853 

1. 39459 

1.  300 

0.60135 

0. 9137L3 

2.  600 

4. 38617 

1. 51426 

1. 250 

0.  50003 

0, 917210 

2.  500 

3.  88115 

1. 68686 

1.  200 

0.  39967 

0.921956 

2.  400 

3.45659 

1. 80551 

1.  175 

0.34957 

0.928299 

2.  300 

3.  08831 

1. 95835 

1.  150 

0.  29959 

0.936595 

2.  200 

2.  75875 

2.  16338 

1.  125 

0.  24961 

0.947424 

2.  100 

2.  46005 

2.45458 

1.  100 

0.  19967 

0.961500 

2.  000 

2.  18679 

2.90595 

1.  075 

0.  14979 

0. 979892 

1.900 

1.92637 

3. 71930 

1.050 

0.09981 

1.00415 

1.  800 

1. 68321 

5. 76345 

1.025 

0.04995 

1.03666 

1.  700 

1. 45191 

Co 

1.000 

0 

Fig.  8.  Any le  of  deflect 


